Kondo universal scaling for a quantum dot coupled to superconducting leads 
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We study competition between the Kondo effect and superconductivity in a single self- assembled 
InAs quantum dot contacted with Al lateral electrodes. Due to Kondo enhancement of Andreev 
reflections the zero-bias anomaly develops sidepeaks, separated by the superconducting gap energy 
A. For ten valleys of different Kondo temperature Tk we tune the gap A with an external magnetic 
field. We find that the zero-bias conductance in each case collapses onto a single curve with A/ZcsTk 
as the only relevant energy scale, providing experimental evidence for universal scaling in this system. 
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The concept of universality evolved from the study 
of phase transitions in the field of statistical mechan- 
ics, where it was found that critical phenomena can be 
grouped in classes. All systems within the same univer- 
sality class show identical behavior near the critical point 
for a certain set of relevant observables, despite the dif- 
ferences in microscopic details of the systems at hand [1]. 
Although formally not a phase transition, the Kondo ef- 
fect has been shown to exhibit universality as well [2]. 
The Kondo effect arises from the anti-ferromagnetic in- 
teraction of localized magnetic impurities with conduc- 
tion electrons. Below a characteristic temperature, the 
Kondo temperature T^, a strongly correlated many-body 
state is formed that screens the impurity and greatly in- 
creases its scattering cross-section [3] . The last ten years 
have seen a renewed interest in this phenomenon with the 
possibility to use the unpaired electron spin on a quan- 
tum dot (QD) as the impurity, a setup that offers un- 
precedented control over experimental parameters [3, [HI . 

Kondo universality manifests itself as the governance of 
the Kondo energy scale ksTx over all low-energy physics 
of the system. Other scales, such as the charging energy 
Uc or lead coupling F, can be an order of magnitude 
larger, yet to describe the response to external param- 
eters we require only Tk and Gq^ the zero temperature 
conductance. Under influence of a finite temperature, 
bias voltage or magnetic field {X = ksT^ eVsD or gjisB) 
the normalized conductance G{X)/Go is reduced mono- 
tonically from unity following a curve that scales only 
with X/ksTx [6]. Once appropriately scaled the behav- 
ior of any microscopic realization of the Kondo model 
can be collapsed onto this curve. The prime example of 
Kondo universality is its temperature dependence, where 
the empirical function 



G{T)/Go= [1/(1 + {T/nf)y 



(1) 



is widely used to fit experimental data. Here 
Tx/V2-'-/^ — 1 , to conform with the definition of the 
Kondo temperature G{Tk) = Go/2, and s = 0.22 for 
spin- 1/2 systems f?\. Once expressed in the dimension- 



less parameters G{T)/Go and T/Tk the behavior of e.ff. 
semiconductor based dots Q, carbon nanotubes [1, 0] 
and single-molecule transistors [lO] could be mapped 
onto curve [H showing the dominance of Tk over the sys- 
tem's microscopic details. 

Both superconductivity and the Kondo effect are 
many-body interactions with a strict inherent spin or- 
dering. In a regular s-wave superconductor the electrons 
condense in singlet Cooper pairs, whereas the Kondo 
ground state consists of a dynamic many-body singlet 
formed between the localized dot electron and the itin- 
erant lead electrons. Their competition is characterized 
by the relative strength A/Zc^Ti^; for A/Zc^T^ <C 1 the 
Kondo effect can survive in the presence of SC leads, in 
the opposite limit A/Zc^^V ^ 1 Kondo is suppressed 
restoring Coulomb blockade (CB). The possibility of 
Andreev scattering at the dot-SC interface results in 
an even richer physics. Through Multiple Andreev 
Reflections (MAR) a Cooper pair of charge 2e can be 
transported across the dot In the Kondo regime a 
resonant channel opens up at the Fermi level, making 
the dot more transparent for the MAR processes. This 
allows the conductance to exceed the unitary limit of 
2e^//i, provided that Tk exceeds A [12]. This surprising 
prediction was confirmed experimentally by Buitelaar et 
al. [9] for a carbon nanotube with Al electrodes. They 
observed a crossover between enhancement and suppres- 
sion of the Kondo conductance at ksTx ^ A, explained 
with a resistively-shunted junction model We work 
in a different regime where MAR is strongly damped 
by repulsive Coulomb interactions, though we will see 
that first order AR is enhanced by the Kondo effect. 
In contrast to reference [o^ the zero-bias conductance is 
suppressed over the entire range of Kondo temperatures 
by the gap A at the Fermi level. 

We shall study transport through a single self- 
assembled InAs quantum dots contacted laterally via the 
nanogap method [14] to Al SC electrodes. The dots are 
formed by MBE deposition of ~ 4 ML InAs on a GaAs 
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FIG. 1: (Color online) (a) Differential conductance dl/dVsD 
as a function of Vg and Vs d , with the leads driven normal by 
a field B = 100 mT. T = 35 mK. (b) dl/dVsD at V^d = mV 
in units of /h. We show KV5 to KV7 with a Tk of 1.04, 1.65 
and 1.18 K respectively, (c) Scanning Electron Micrograph of 
a similar device as used in these studies, with 100 nm scale 
bar. An InAs dot bridges the source and drain electrodes, (d) 
Temperature dependence of the conductance at KV5, with a 
best fit to scaling curve [H 



substrate, at this coverage many individual dots have co- 
alesced into larger InAs islands with a diameter of ~ 80 
nm. The source and drain contacts, separated by a 25 nm 
gap, are defined by e-beam lithography. E-beam deposi- 
tion of Ti (5 nm) and Al (100 nm) is preceded by a 5 s 
BHF etching step to de-oxidize the surface. Roughly 5% 
of the gaps is bridged by a dot, yielding a working device 
(fig. [It). Each Al electrode is doubly contacted, allowing 
4- terminal measurements. A degenerately doped layer 
300 nm below the surface acts as an electrostatic back- 
gate. Using a 4-terminal Al test strip deposited simulta- 
neously we can characterize the Al leads to have a critical 
temperature Tc — 0.9 K and critical field Be = 95 mT 
in perpendicular orientation. More details of the fabri- 
cation process will be given elsewhere [15]. 

The measurements presented were done on a device 
with a RT junction resistance of 100 kQ. First we ap- 
ply a 100 mT field to suppress superconductivity in the 
electrodes. Figure [1^ shows the charge stability dia- 
gram (Coulomb diamond) of the device, where we control 
the number of electrons on the dot with the voltage Vg 
applied to the back gate. Figure [TJd shows the linear 
response at zero bias. We find a dot charging energy 
Uc ~ 2.5 meV and an electronic level spacing 5E ~ 2.5 



meV. Counting the Coulomb oscillations we obtain an 
electron occupancy of = 12 for the leftmost valley at 
Vg = 3.0 V. Around = 3.64 V we reproducibly suffer 
from an electrostatic switch in the environment, which 
we attribute to the charging of a nearby InAs dot. 

The stability diagram clearly exhibits the even-odd 
parity associated with spin-1/2 Kondo, showing an en- 
hanced conductance ridge around zero bias for odd occu- 
pancy. For the various Kondo valleys the value of Tk is 
different each time, as the electronic orbitals couple dif- 
ferently to the reservoirs. In addition to this variation, 
thermal cycling to room temperature also alters the dot 
parameters including Tk- Divided over three cooldowns 
of the device we were able to study 10 different values of 
Tx, ranging from 0.7 to 2 K. We shall refer to them as 
Kondo Valley (KV) 1-10. For KV5 the zero-bias conduc- 
tance is shown as a function of temperature in figure [T]l, 
the line shows a best fit to the empirical curve [TJ The 
fitting gives Tk = 0.98 K, which is close to the value 
of 1.04 K found from half the width of the Kondo reso- 
nance out of equilibrium. From now on we shall use the 
latter method to determine Tk- The maximum normal 
state conductance Gq was always smaller than 2e^//i, 
indicating an asymmetric lead coupling. 

Sweeping back to zero magnetic field we bring the 
electrodes in the SC state. Figure [2^ shows the zero field 
Coulomb diamond. For even occupancy two parallel 
lines at source-drain (SD) bias \Vsd\ = 2A/e mark the 
onset of direct quasiparticle tunneling. Their separation 
of 4 A allows us to determine A = 0.14 meV, in good 
agreement with the critical temperature of the leads 
via the BCS relation A(T = 0) = Ll^ksTc [le^. An 
accurate scan of the gap structure at even electron 
number is shown in Fig[2)3, the peaks at \Vsd\ = 2A/e 
reflect the singularity in the density of states (DOS) 
of the leads. The MAR found in other QD devices 
connected to SC-leads 17, lisl is here absent due to 



strong on-site Coulomb interaction between the carriers 
(Uc ^ 17 A) [ll|. Only the first order Andreev reflection 
at I^/sdI = A/e is faintly visible as a shoulder to the 2 A 
feature. 

We now turn to the behavior in the Kondo regime. 
The transition from the N- to SC-state with the magnetic 
field is shown in Fig [2]3. With the leads in the N-state 
the field splits the Kondo resonance; we observe a peak in 
the differential conductance when Vsd equals the Zeeman 
energy A linear fit gives a g-factor \g\ = 6.6 ±0.3. 
Below Be we surprisingly observe three resonances in 
the differential conductance. The two non-equilibrium 
ones occur at \Vsd\ = A/e, reflecting the SC DOS of 
the leads. Their position follows the field dependence of 
A, being maximum at zero field. This conspicuous posi- 
tioning corresponds to the first order Andreev reflection 
(AR), the physical process where an incoming electron 
scatters at the N-SC interface, injecting a Cooper pair 
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FIG. 2: (Color online) (a) dl/dVsD as a function of Vg and 
VsD, with the leads in the SC-state {B = OT). (b) Scan of the 
gap structure in the CB regime, as indicated by the square in 
upper figure. Conductance is suppressed for Vsd < 2 A. (c) 
Color-scale plot of dl/dVsD, showing the field dependence of 
the Kondo resonance (KV2, Tk = 1.23 K). Horizontal lines 
denote the SC-N transition in the leads, (d) The Kondo reso- 
nance out of equilibrium for three different Kondo valleys (KV 
10,2,4 from upper to lower) with the leads in the SC state. 
The curves are normalized by their N- state conductance Gq . 
Offset for clarity, (e) Schematic depiction of the Kondo en- 
hanced Andreev Reflection, with a split Kondo resonance in 
the dot DOS. (f) Zero-bias Kondo conductance (KV2) as a 
function of magnetic field, with N-state leads. 



in the SC-region and a backreflected hole in the N-region 
[l6] . Trough an interaction with the Kondo effect the AR 
is enhanced [13., 20,]; the mechanism is depicted schemat- 
ically in Fig [2^. In the middle of the valleys charge fluc- 
tuations are inhibited and the electrons traverse the dot 
through cotunneling processes. In the CB regime the 
lowest available state is ~ {Uc + SE)/2 away, but, by 
contrast, in the Kondo regime the many-body resonance 
gives rise to a finite DOS at the Fermi energy, enhancing 
the conductance. At finite bias the Kondo DOS splits 
with a maximum at each lead chemical potential, for 
Vsd = A/e this maximum aligns exactly with the band 
edge of the opposite reservoir. The electrons that con- 
tribute to the AR tunnel at this energy, and in this way 
the first order AR is strongly enhanced by the Kondo 
effect. 

In Fig [2]i we compare the normalized conductance 
G^^ /Gq as a function of Vsd for 3 different valleys. 



Surprisingly the height G^^\a of the satellites appears 
to be insensitive to A/ksTK^ and was between 0.65 and 
0.9Gq for all 10 values of Tk studied. A very recent 
work by Sand-Jespersen [21[ reports on similar Kondo 
enhanced AR in dots of small Tk- Why this process can 
persist where A ^ ^bTk is uncertain. We speculate 
that the exact alingment with the band edge allows an 
out-of-equilibrium recovery of the Kondo state similar to 
the case of an applied magnetic field. Note that the di- 
rect tunneling onset at \Vsd \ = 2A/e is not enhanced by 
Kondo, and its features are dwarfed by the much larger 
A-peaks. 

In contrast to the A-satellites the zero-bias conduc- 
tance is very sensitive to the relative strength A/Zc^Tx, 
or the ability of the Kondo effect to break up lead 
spin-pairing. For A I^bTk (upper graph) the zero-bias 
conductance is close to its N-state value, whereas for 
higher /S./kBTK the SC gap suppresses the conductance 
increasingly. Note that the resonance around Vsd = ^ 
cannot be interpreted as a supercurrent, as was done 
by Grove- Rasmussen et al. [22]. A low critical current 
Ic ~ 100 pA is to be expected as the Cooper pairs of 
charge 2e cannot tunnel directly due to the high Uc- 
The lead noise temperature of ~ 100 mK is roughly 10 
times higher than the Josephson energy Ej = hlc/2e, 
destroying SC phase coherence across the junction. 
Indeed no dissipationless fiow was detected in a current 
biased 4-terminal measurement. 

We analyze the 10 separate Kondo valleys with each a 
different Tk- On top of that we use the magnetic field to 
alter the SC gap, providing us a knob to experimentally 
tune A/ksTK' By comparison, using the Vq dependence 
of Tk as a knob gives a much smaller range of variation. 
In Fig [3] we have plotted G^^ /Gq at diflFerent fields as 
a function of A/kBTK- When plotted in this fashion 
all data points collapse onto a single curve, proving that 
A/kBTK is the only relevant scaling parameter for the 
system. For comparison the inset shows the same data 
set as a function of just Tk- The scaling curves for other 
external parameters, such as temperature [23^], bias [IJ] 
or irradiation [iH (X = /c^T, eVsD and huo)^ are charac- 
terized by a logarithmic conductance fall-off in the weak 
coupling regime {X kBTK) where a perturbative ap- 
proach is valid. Although the data range is too limited to 
unambiguously claim any functional form, the observed 
fall-off seems to agree with the expected logarithmic de- 
pendence. Fitting the data with an empirical form such 
as function [T] gives a poor correspondence. However, the 
kBT and A scaling curves need not be the same as the 
underlying mechanism is fundamentally different. Fur- 
thermore we believe the shape of the curve to be slightly 
distorted by the magnetic field, as discussed later on. In 
contrast to the results of Buitelaar we always find that 
the gap A suppresses the zero bias conductance, also in 
the regime kBTK > A . The main difference is the strong 
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FIG. 3: (Color online) Semi-log plot of G^^ /G^ for KVl-10 
at various fields. Plotted against A/ksTK the data collapses 
on a single curve. The field dependent value of A was 135, 
140, 120, 100 and 60 /xeV for 0, 20, 40, 60 and 80 mT respec- 
tively. The measurement uncertainty is given by the bars. 
Inset shows the same data plotted against Tk only. 



Coulomb repulsion in our device that damps the MAR, 
the latter being necessary to raise G^^ above Gq . Un- 
fortunately the field behavior is not discussed in reference 
[qI, making a direct comparison difficult. 

To shed some light on the origin of the observed scaling 
we want to point out a similarity with the effect of other 
perturbations. Under equilibrium conditions the spin-flip 
cotunneling processes that create the Kondo singlet con- 
serve energy. At finite Vsd or ksT there is phase space 
available to inelastic processes J^ading to decoherence of 
the many-body Kondo state [3, [25[. In the case of SC 
leads a spin-flip with the leads requires the breaking of a 
Cooper pair, creating two excitations in the conduction 
band. This cost of 2A requires energy exchange with 
the environment, leading to a suppressed conductance 
through decoherence of the Kondo singlet. 

Finally we want to address the non-negligible Zeeman 
splitting caused by the field that tunes A. At 80 mT it 
equals Ez =30 /ieV, which should be compared to the 
Kondo energy ksTx that ranges from 70 to 170 /ieV. 
The suppression of the conductance by the splitting 
depends on both \B\ and T^, and consequently we can 
expect the shape of the observed scaling function to 
be slightly distorted. Furthermore we used G^ at 100 
mT, rather than the 'true' Go, for normalization. In 
figure [2f we have plotted the equilibrium conductance 
of KV2 as a function of magnetic field, with SC-leads 
regime omitted. In the absence of a proper scaling curve 
for the 5-dependence we have fitted a parabola, and 
estimate that Go is 12% higher than Gq. Because the 
field reduces both G^^ and Gq the 12% represents an 



upper bound on the error. This is still smaller than the 
measurement uncertainty as indicated by the bars; we 
therefore believe the influence of the unwanted Zeeman 
splitting does not falsify our claim of Kondo universal 
scaling. 
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